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We study response and velocity autocorrelation functions for a tagged particle in a shear driven 
suspension governed by underdamped stochastic dynamics. We follow the idea of an effective con- 
finement in dense suspensions and exploit a time-scale separation between particle reorganization 
and vibrational motion. This allows us to approximately derive the fluctuation-dissipation theorem 
in a "hybrid" form involving the kinetic temperature as an effective temperature and an additive 
correction term. We show numerically that even in a moderately dense suspension the latter is 
negligible. We discuss similarities and differences with a simple toy model, a single trapped particle 
in shear flow. 

PACS numbers: 82.70.-y, 05.40.-a 



I. INTRODUCTION 

Equilibrium statistical mechanics describes the connec- 
tion of a few macroscopic intensive quantities, e.g. pres- 
sure and temperature, with the microscopic properties of 
the many particles constituting the system [T] . Although 
there is no equivalent formalism for systems driven out of 
thermal equilibrium studies are often motivated by the 
quest for simple governing principles such as an effec- 
tive temperature f?. "iP. Of particular interest are small 
mesoscopic systems such as colloidal particles, nanopar- 
ticles in solution, or biological systems, all of which are 
dominated by fluctuations. 

For systems only slightly perturbed from equilibrium 
- into what is called the linear response regime - the 
fluctuation-dissipation theorem (FDT) relates the re- 
sponse with equilibrium correlations through the tem- 
perature 2]. This unique temperature corresponds to 
what we would measure with a thermometer, and is in- 
dependent of the observables entering the fluctuation- 
dissipation theorem. Examples for these observables are 
the velocity related to the diffusion coefficient, or stress 
fluctuations determining the viscosity. The practical im- 
portance of the FDT both for experiments and simula- 
tions thus stems from the possibility to extract transport 
properties from stationary fluctuations. In a more re- 
cent application the FDT has been used to predict the 
self-assembly of model systems [5, 3^ . 

Even for systems driven into a non-equilibrium steady 
state (NESS) far beyond the linear response regime we 
can still define a linear response function [THlOj. For 
the broad class of driven systems with Markovian dy- 
namics and embedded in a fluid at well defined temper- 
ature it has been demonstrated recently that then the 
FDT is most convincingly interpreted in terms of an ex- 
cess correlation |11H17| . The single temperature that 
enters these generalized FDTs is that of the fluid, and 
no approximations are involved. Nevertheless, for the 
purpose of a simple description we might still be inter- 



ested in defining an approximate temperature through 
the fluctuation-dissipation ratio (FDR), i.e., the ratio 
of correlation to response function. This strategy has 
originally been proposed in the context of aging mean- 
field spin systems |3] , and subsequently been applied 
to many different systems [rSlEU] . In particular, in a 
sheared colloidal suspension or fluid the Einstein relation 
between the self-diffusion coefficient of a tagged particle 
and its mobility is broken and can be used to define an 
effective temperature pTM25j . 

In this paper we study a many-body system governed 
by underdamped stochastic dynamics and driven into a 
NESS through linear shear fiow. Such a system could 
model a fiuid with every particle coupled to a stochas- 
tic thermostat, colloidal or nano-suspensions, or dusty 
plasmas We follow a tagged particle, i.e., a ran- 

domly chosen particle out of many identical interacting 
particles. Motivated by the physical picture of an ef- 
fective confinement in dense systems we also consider a 
single trapped colloidal particle in shear flow [37] as a 
toy model, see Fig. [l] We discuss the FDT in a "hybrid" 
form: we relate response and correlations through the 
kinetic temperature in the spirit of an effective temper- 
ature, but with an additive correction term still present. 




FIG. i: (Color online) Systems studied: (A) Toy model with 
a single particle bound to the origin and (B) tagged particle 
in a dense suspension. Both systems are driven through linear 
shear flow. 
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For the trapped particle expressions can be obtained an- 
alytically. For the tagged particle we derive a similar 
FDT exploiting a time-scale separation due to the eflFec- 
tive confinement. In both cases we show that the correc- 
tion term indeed becomes negligible for strong confine- 
ment. 



II. FLUCTUATION-DISSIPATION RATIO 

We study the response of a single particle with mass 
M moving in a viscous liquid at temperature T. After 
applying a small external force f directly to the particle 
its mean velocity evolves as 



(v(0) = / dt' R(< 
Jo 



t'){{t') + 0{f 



(1) 



where R = {Rij ) is the response matrix with components 



Rij{t — t') 



(2) 



f=0 



with t ^ t' due to causality. The brackets (• • •) denote the 
thermal average in the perturbed system. Throughout 
this paper we employ dimensionless units and measure 
length in units of the particle diameter a and energy in 
units of k^T . Time is measured in units of tq = / Dq, 
which quantifies the time it takes for a free particle with 
diffusion coefficient Dq to diffuse a distance equal to its 
diameter. The reduced mass m = {MDQ/k-QT)/Ta is the 
ratio of the momentum relaxation time to the diffusive 
time scale tq. 

We describe the stochastic particle motion through the 
coupled equations r = v and 



-V[/ + f - [v - u(r)] + ^, 



(3) 



where r and v are the particle position and velocity, re- 
spectively. Besides the conservative forces arising from 
the potential U we can perturb the particle by a direct 
force f . The noise ^ modeling the interactions of the 
particle with solvent molecules has zero mean and corre- 
lations 



(IWI^(0)=2 1<5(t-0- 



(4) 



The external shear flow enters through the term u(r) — 
jyex, i-e., the flow points in x-direction and increases 
its amplitude linearly with the y-coordinate. Here, 7 is 
the strain rate which in our units is equal to the Peclet 
number. 

Eq. ^ measures the linear response of the system. If 
the unperturbed system (f = 0) is in thermal equilibrium 
the FDT 



(5) 



relates this response to the velocity auto correlation func- 
tion (VACF) C{t), where the subscript indicates that 



these correlations are to be measured in the unperturbed 
system. 

For a computationally convenient representation of the 
response function Eq. ([2| valid both in and out of equi- 
librium consider the path weight of the noise 



(6) 



The stochastic velocity of the particle is a result of pre- 
vious collisions with solvent molecules, v(t) = v[f;^(r)]. 
Since a force perturbation is equivalent to perturbing the 
noise we can write 



R,,{t-t') 



(7) 



for the response. A functional integration by parts then 
leads to [ni [55] 



R(^-0 = ^(v(tK^(i')>o. 



(8) 



Hence, even if the system is driven into a NESS the re- 
sponse can be measured through a steady-state correla- 
tion function (see also Refs. [291 [30]). However, the FDT 
in the form Eq. ([5| no longer holds and the dimensionless 
fluctuation-dissipation ratio (FDR) is defined as 



Rii(t) 



(9) 



for the diagonal components. 

To calculate Xi{0) we perform a short-time expansion 
of Eq. ^ leading to (v) w AtR(0)f. From the equa- 
tion of motion ([s]) we obtain m(v)o ~ m{v) / At — f. 
We have exploited that the average force on the parti- 
cle right before the perturbation vanishes, (VJ7)o = 0, 
which is quite obvious for isotropic systems but due to 
the inversion symmetry about the origin it holds also in 
the presence of simple shear flow. Hence, in our dimen- 
sionless units we obtain Rij{0) = Sij/m and finally 



(10) 



The right hand side is the kinetic temperature as mea- 
sured through the velocity fluctuations. In the following 
we study for two systems whether, and under which con- 
ditions, Eq. (10) extends to times t > 0, i.e., whether 



III. TRAPPED PARTICLE 

The toy model we investigate first is a single particle 
trapped in the harmonic potential 

C/(r) = ^-kr' (11) 

with strength k, where r is the displacement from the 
origin and r — |r|. Due to the linearity of the restoring 
force the z-component in Eq. (|3| decouples and remains 
in equilibrium, Xz{t) — 9z — \. Therefore, in this section 
we only consider the motion in the xy-plane. 
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A. Analytical results 



Due to the quadratic potential Eq. (11) the equations 
of motion comprise a Unear system of first order differ- 
ential equations. Hence, we can solve it for the velocity 

v(i) = G^'-(t)ro + G™(i)vo + - / di' G™(i-i')C(i') 

™ Jo 

(12) 

given the initial displacement Tq and initial velocity vq. 
The explicit expressions for the Green's functions G(t) 
are given in the appendix |X] Both the VACF and the 
response function are easily calculated from the solution 
Eq. For the VACF we find 



C{t) = (v(t)vJ>o = G-(t)(rovT)o + G™(t)(vovJ)o, 

(13) 

while the response function is trivially related to the 
Green's function through 

R{t - t') = i(vWI^(i'))o - ^G™(< - t') (14) 

using the noise correlations Eq. Q. From the steady 
state distribution Eq. (All we obtain the moments 



(i-qVo )o = 



1 

2k 



(vovo)o 



1 ^ 1 

m 2k 



-7 
7 

72 




(15) 
(16) 



The kinetic temperatures Eq. ( 10 ) perpendicular to the 
shear flow are 6y = t 



1, whereas 



1 + a^r > 1, 



m 
2k' 



(17) 



As a consequence of the linear forces and the symmetry 
of the shear flow the excess compared to equilibrium is 
proportional to 7^. 



B. The FDT 

Using the explicit expressions for the moments we 
see that with Gy'^{t) = the correlation function for 
the y-component reads Cyy{t) = G™{Vy)a and therefore 
Xy{t) = 9y = 1 at all times. On the other hand, in the 
direction of the shear flow we obtain 



C^t) = 9,,R,,{t) + {yv,)oGllit), 



(18) 



i.e., the velocity correlations are expressed through the 
response times the kinetic temperature plus a correction 
term. Separating the dependence on strain rate the cor- 
rection term can be rewritten as 



(19) 



where max|/j:(t)| — 1 and f3x captures the magnitude 
of the correction term. In Fig. [2] we plot scaled re- 
sponse and correlation functions for three representative 
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FIG. 2: (Color online) Single particle moving in a har- 
monic trap: (A-C) scaled velocity auto-correlation function 
Cxx{t)/9x and response function Rxx(t) vs. time t for strain 
rate 7 = 2 and different masses m and trap strengths k. 
(B,C) For large k, the correction in Eq. ( 18 1 becomes neg- 
ligible and both curves lie on top of each other. (D) The 
magnitude of the correction term as a function of k for 
two different masses. (E) Sketch of the different regimes of 
the FDT for m — 0.6. The dashed lines 7^ and 1/m limit 
the region where the FDR is approximatly time-independent, 
Xx ~ Ox- Below the solid line 7717^ the kinetic temperature 
is much larger than unity. While for the chosen m there is 
a gap (shaded area), with increasing m both regimes can be 
realized in the vicinity of the solid line. 



values of the parameters m and k. As demonstrated in 
Fig. [2JI) increasing the trap strength k strongly decreases 
The behavior of the FDT is sketched in Fig. [2p, 
where we compare the magnitude of the correction term 
to the kinetic temperature. From Eq. (17 1 we see that 
for k <^ TO7^ we have 9^ ^ I (see also Fig. [2^ 



For the 

relevant case k > 1/m we find from the explicit expres- 
sions that for the FDR to become approximately time- 
independent k ^ max{7^,m7^} must hold (Fig. |2]3). 
While for small masses m < 1 there is a gap (this case 
is sketched in Fig. [2J1;), for sufficiently large m > 1 these 
two regimes come close. In Fig. [213, we demonstrate that 
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there is indeed a regime of intermediate trap strength 
k ~ mj'^ with an increased effective temperature where 
nevertheless X^it) ~ 9^ holds to a very good degree. 



IV. TAGGED PARTICLE IN A SUSPENSION 

The main system we now study is a suspension com- 
posed of N particles in which we tag and follow a single 
particle, say fc = 1, with position ri and velocity Vi. The 
particles interact through a pair potential u{r) with to- 
tal potential energy U = J2i<j'^{\^i ~ Advection 
through the shear flow leads to the well-known Taylor 
dispersion [3T] . Instead of the absolute velocity Vi it will 
be more convenient to use the relative velocity 



u(ri) 



(20) 



with respect to the flow as the observable entering the 
response R{t) and correlation function C{t). 



A. Time-scale separation 

In the following we assume a time scale separation be- 
tween the motion of the potential energy minimum (or in- 
herent state position ^32^ ) r^, and the vibrational motion 
of the tagged particle around Vc- The physical picture is 
that particles vibrate in a "cage" of surrounding parti- 
cles and that local reorganization takes much longer than 
the vibrational motion. Linearizing the force exerted by 
neighboring particles on the tagged particle leads to 



-Vi[/«-k[ri(i)-re(t)], 



(21) 



We solve the resulting equations of motion leading to 



the same formal result Eq. ( 13 ) for the correlation func- 
tion and Eq. (14) the response function. In principle 



the Green's function can be calculated but wc will not 
need its explicit form here. To see that the flrst term in 
Eq. ( 13 ) vanishes consider the projected probability 



V's(ri,v) = y"dr2---drAr y"dv2' 



■dvA, V^({rJ) (22) 



of the tagged particle, where ips is the full stationary dis- 
tribution. For a homogeneous, translationally invariant 
system cannot depend on the position ri and therefore 
(v)o = vanishes. Hence, we obtain 



(23) 



as our central result. The correction term differs from 



Eq. ( 18 1 and now couples to the off-diagonal elements 



of the response function instead of the off-diagonal com- 
ponent of the Green's function connecting velocity and 
position. 



B. Langevin dynamics simulations 

We perform Langevin dynamics simulations to study 
Eq. (23) for a specific system. The N — 1728 particles 



are enclosed in a cubic simulation box with edge length L. 
The particles interact through the Yukawa (or screened 
Coulomb) pair potential 



u{r) = 



(r^l) 
(r < 1), 



(24) 



where e is the interaction energy at contact and the 
screening length determined by the composition of the 
surrounding solvent. We choose e = 8.0 and k = 5.0 in 
order to obtain a broad range of densities for which the 
liquid phase is stable [33]. For the shear flow we em- 
ploy Lees-Edwards boundary conditions, which enforce 
a linear velocity profile in the suspension |34j . We in- 
tegrate the equations of motion by a stochastic velocity 
Verlet algorithm |3S] with a time step 5 x lO"**. Since 
the hard-core repulsion cannot be implemented in the 
interaction potential we employ a simple algorithm that 
detects collisions and computes the appropriate positions 
and velocities after the impact according to momentum 
and energy conservation (see Refs. '3S', 'ST' and references 
therein). The NESSs are prepared by initializing the par- 
ticle positions on a regular lattice at low density. Then 
we equilibrate the system and slowly increase the density 
by decreasing the box size. After this equilibrium system 
is constructed we slowly ramp up the strain rate until 
the target value 7 is reached. We simulate another 1000 
time steps to relax the system into the steady state. This 
procedure is repeated separately for independent runs. 

We study suspensions at different volume fractions 
(j) = TrN/{6L^) and strain rates 7. To determine the 
response and correlation functions we simulate the mo- 
tion of the system in the NESS and record the velocity 
trajectories of 200 randomly chosen particles in four in- 
dependent runs. From this data we can easily evaluate 
the VACF. To compute the response to a small force we 
again employ Eq. ([s]), which is still valid for the many 
particle system. This enables us to obtain the response 
function from steady state trajectories without the need 
to explicitly perturb the system. Therefore, in addition 
to the velocity, we also record the stochastic forces acting 
on the tagged particles, which are directly accessible in 
a numerical simulation. 



C. Numerical results 

The pair distribution function quantifies the probabil- 
ity to find another particle at a displacement r. It is dis- 
torted from its isotropic equilibrium shape, see Fig. |3]4.. 
This function visualizes the average environment of the 
tagged particle. Under the influence of the shear flow 
there is a higher number of particles in the compressional 
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FIG. 3: (Color online) Tagged particle with m = 1 in a sheared suspension at volume fractions (j> — 0.1 (top row) and tj) = 0.4 
(bottom row): (A) Pair distribution function in the xy-plane for strain rate 7=1. (B) Scaled velocity auto-correlation functions 
Cxx{t)/9x and response functions R^xit) for different strain rates: from bottom to top 7 = 0,0.2,0.4,0.6,0.8, 1. For visibility 
curves are shifted. (C) Off-diagonal response functions Rxy{t) (solid lines) and Ryxit) (dashed lines) for the same strain rates 
and offsets. (D) The magnitude of the correction term \6xyRij\, i.e., the product of Eq. (25 1 with the maximum i?*^ of the 
off-diagonal response function, as a function of strain rate. The dashed lines show the quadratic fits. 



zones and a depletion in the extensional zones p7] . More- 
over, the peak corresponding to the first nearest-neighbor 
shell becomes more pronounced. 

In Fig.[3j3, the response functions Rxx{t) together with 
the scaled VACF Cxx{t)/9x is plotted for volume frac- 
tions = 0.1 and (j) = 0.4, and for a range of strain rates. 
Even for times i > these functions coincide, which im- 
plies that the additive correction term in Eq. ( 23 1 van- 



ishes. The off-diagonal response functions are plotted in 
Fig. [3p. For comparison response and correlation func- 
tions are plotted in Fig. |4] for a dilute suspension at vol- 
ume fraction (j) = 0.01. Here we observe a clear difference 
between correlations and response. To understand the 
observed behavior one should bear in mind that we con- 
sider the velocities relative to the local flow. At low den- 




FIG. 4: (Color online) Dilute suspension at = 0.01 and 
strain rate 7=1 (left) and 7 = 2 (right). In comparison 
to Fig. [3] a deviation between response Rxx (t) and velocity 
autocorrelation function Cxx{t) is observed. 



sities collisions are rare and the particles adapt smoothly 
to the flow of the solvent. Therefore, the largest devia- 
tions from the flow proflle arise when a particle diffuses in 
y-direction entering a region of faster or slower flow in x- 
direction, to which it needs to adapt. For higher densities 
collisions become much more frequent. These collisions 
prevent the particles from adapting to the solvent flow. 
Additionally they distribute the momentum, transfered 
to the particles by the shear flow, more or less randomly 
in the three spacial directions. The result is that under 
shear flow the diagonal velocity moments grow with in- 
creasing density, and become more and more similar in 
size. 

The two components {vxVz)o ~ and {vyVz)o ~ of 
the velocity correlation matrix are very small (see also 
Fig. [t]) . Defining the off-diagonal "temperature" 



•ixy = ■m{VxVy)o 



(25) 



we see that the dominant contribution to the correction 
term in Eq. (23 1 is OxyRxy(t) for the x component, and 



-' xyRy x \ 



r{t) for the y component. In analogy to the trapped 
particle [Eq. ( 19 1] we separate the strain rate dependence 



of the correction terms, 



^xyRxy (^) 



YPxIx{t), 9xyRyx{t) 



^i^(iyly{t), (26) 



with again max \ Ii{t)\ = 1 and coefficients Pi. In Fig.[3jl), 
we plot \9xyR*j \ ~ 7^/3i, where R*j is the maximal value 
of the off-diagonal component of the response matrix. For 
(/) = 0.1 we observe the predicted quadratic dependence 
on the strain rate 7. For (j) = 0.4 the quadratic prcdicition 
holds for 7 ^ 0.8, while for larger strain rates higher order 
terms in 7 become important. 
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FIG. 5: (Color online) The kinetic temperatures 9i (symbols) 
and their quadratic fits (lines) plotted vs. strain rate 7 for 
m = 1. Shown are the three directions for volume fractions 
(A) ^ = 0.1 and (B) (j) = 0.4. In (C) the coefficients a; and 
Px X lO'^ are plotted as function of the volume fraction (f). 
(D) Relative magnitude of the correction term \OxvRxy\/(^x x 



10 as function of strain rate and density, cf. Fig. |2p. Colors 
are linearly interpolated. 



In Fig. [5] we plot the kinetic temperatures as a func- 
tion of strain rate. While for = 0.1 a clear difference 
between motion parallel to the flow (0^) and motion per- 
pendicular to the flow {9y ^ Oz) can be seen, this dis- 
tinction is diminished at higher densities. Moreover, all 
kinetic temperatures can be well fitted by the quadratic 
function Eq. (17) with coefficients for the three di- 
rections. 



The increase of velocity fluctuations can be 
explained by forced collisions due to the flow gradient, 
an effect that is more pronouned at higher densities and 
higher strain rates. 

In Fig. [5p, the coefficients ai and /S^ are shown for 
the different densities. The coefficient f3x decreases for 
larger densities but then turns up again at (jj — 0.4. The 
reason for this non-monotonic behavior is that there are 
two effects determining the shape of the a;?/-component 
of the response function, see Fig. [6] One dominates for 
lower, one for higher densities. At low densities forc- 
ing the particle upwards in y-direction moves it into a 
region of faster flow. While the velocity relaxes the rel- 
ative velocity with respect to the shear flow is negative. 
At higher densities this effect is weaker because of more 
frequent particle collisions, and thus large excursions in 
y-direction are rare. Hence, the velocity differences due 
to the motion are smaller and so is the response caused 
by this effect. In addition another effect of the collisions 
becomes more signiflcant. Pulling the particle upwards in 
y-directions makes collisions with particles from the left 
more likely than with particles from the right. This leads 
to an average acceleration to the right which counteracts 
the first effect. In the intermediate density regime these 
two effects almost cancel and lead to a small fS^- 

Finally, in Fig. Ist) the magnitude of the correction 
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FIG. 6: (Color online) Sketch of the time evolution after 
a force perturbation in j/-direction (arrow) determining the 
shape of the Rxyit) response function: (A) At low densities 
the tagged particle is slower than the surrounding flow field 
due to inertia. (B) At high densities collisions with neighbor- 
ing particles are more likely, pushing the tagged particle in 
the direction of the flow. 



term with respect to the response function, \QxyR*xy\/()x, 
is plotted, see also Fig. [2jil. Note that over the whole 
parameter range studied here this value is lower than 
0.02, i.e., Xx ~ Ox holds to a very good degree. However, 
we see that for large strain rates and low density this 
ratio grows by an order of magnitude. Due to the effects 
described above for the highest density the range of strain 
rates for which the correction term is negligible shrinks 
again. 



D. Overdamped limit 

Our results have been derived for systems with under- 
damped stochastic dynamics. Of greater practical im- 
portance in colloidal suspensions is the overdamped limit 
corresponding to neglecting inertia, m — )■ 0. In Fig. [7] the 
dependence of the coefficients on the reduced mass m is 
plotted. As expected for smaller m the kinetic temper- 
atures approach unity, 9i ^ 1. The additive correction 
term is no longer negligible as Px and Py grow. We thus 
recover the fluctuation-dissipation theorem as described 
in the introduction, in which the bath temperature en- 
ters and the equilibrium form of the FDT is completed 
by an excess correlation function. 

Inserting the overdamped equation of motion for the 
tagged particle into Eq. ([8| one arrives at 



C,,{t)^2R,,{t) + {Fl'\t)F'.'\Q)), 



(27) 



for t > and components i,j ~ y,z perpendicular to 
the shear flow [21]. Here, F*^^' — — ViC/ is the force on 
the tagged particle exerted by its neighboring particles. 
The time-integrated version of Eq. ( 27 ) has been obtained 



previously within mode-coupling calculations 22 , 23 , 38] . 
For systems in which the force-force correlations can be 
neglected Eq. ([27]) predicts a universal FDR Xi = 2. 
It has been argued that such an approximation might be 
justified in dense suspensions close to the glass transition 
under sufficiently large shear. 
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FIG. 7: (Color online) The coefficients Pi (solid lines, left axis) 
and the kinetic temperature dx (dashed line, right axis) as a 
function of the reduced mass m for volume fraction (f> — 0.4 
and strain rate j = 2. 



V. SUMMARY 

We have studied the relation between velocity autocor- 
relation and response function of a tagged particle mov- 
ing in a suspension that is driven into a non-equilibrium 
steady state through simple shear flow. Under the as- 
sumption of a time-scale separation between vibrational 
motion and local reorganization the tagged particle ef- 
fectively behaves like a trapped particle. The diagonal 
components of the tagged particle's velocity autocorrela- 
tion function are then given by 



with expansion of the kinetic temperature 



.(7) 



1 



■ a»7 



(28) 



(29) 



Here, Rii{t;'y) are the response functions Eq. (|2|, Ii{t;'^) 
are functions of order unity, and ai and Pi are coefficients 
independent of the strain rate 7. While these expres- 
sions are exact for the trapped single particle our numer- 
ical results show that they hold approximately to a very 
good degree for a tagged particle moving in an interact- 
ing colloidal suspension. One might of course anticipate 
the quadratic dependence on strain rate from symme- 
try arguments close to equilibrium. We have shown here 
that these expressions follow from a time-scale separation 
caused by an effective confinement. 

The effect of the shear flow is to break symmetry and 
to couple the particle velocity to earlier perturbations 
perpendicular to its velocity. This leads to an additive 
correction term that grows with 7^. Since /3i <C 1 this 
correction is negligible up to dimensionless strain rates 
7 ~ (max{/3i})~^/^, which can be far from equilibrium. 
Specifically, here we have studied strain rates in the range 
7^1 and found excellent agreement between correlation 
and response functions, see Fig. [Sj However, already at 
7 > 2 the kinetic temperatures start to divert from the 
quadratic law, indicating the importance of higher or- 
der terms. Increasing the density the tagged particle 



interacts more strongly with its surrounding particles. 
The kinetic temperature increases due to more frequent 
collisions with neighboring particles in conjunction with 
transport due to the flow. This is in contrast to the 
trapped particle where tightening the trap reduces fluctu- 
ations and therefore the kinetic temperature approaches 
the bath temperature. 

An intriguing perspective is to apply our results to su- 
percooled (or supersaturated) conditions. Since dynam- 
ics slows down dramatically and the time-scale separation 
between vibrations and long-lived particle displacements 
becomes even more pronounced we expect that our re- 
sults extend into the supercooled regime. While we have 
employed stochastic dynamics one might speculate that 
our results also hold in systems governed by determinis- 
tic dynamics such as the SLLOD equations of motion [39] 
as employed in Ref. [21]. Future work will also address 
the influence of hydrodynamic interactions and flow gen- 
erated through boundaries. 
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Appendix A: Green's function for a particle in a 
harmonic trap in shear flow 



The equations of motion ([3]) for the potential Eq. ( 11 ) 
are most conveniently written as x = Ax -I- (0, ^)"^ for a 
vector x = (x, y, Vx, Vyj^ with 



A = 









m 














m 


-k 


7 


-1 








~k 





-1 



The Green's function is 



G{t) = e^* = 



8 



We need the explicit expressions for the foUowing two 
matrices: 



1 



sinojil + 'ygyrit)la 



_ {Am?Lj'^ — l)sina;t + (4m^w^ + l)wicosa;i 



cos ujt 



sincji 
2muj 



J-xy 



9vv{t) 



_ {imuP't — l)s,m. ujt + ijjt COS Lot 



where 



UJ = 



\/Akm — 1 
2m ' 



>-xy 



1 




The diagonal components are independent of the strain 
rate. To determine (ix as defined in Eq. ( 19 ) we calculate 
numerically the maximum of G'"Jy{t). 

For the sake of completeness, the other two matrices 
are 



- p-t/2" 



smujt , 

COSut + 1 + -fgrr{t)lxy 



9rr{t) 



_ {2muP't + 1) sin ujt — ujt cos tut 



t/2m 



grv (t) 

The stationary distribution 
1 



Am?u!^ 
sinwt 



1 + igrv{t)ix 



sin cot — Lot cos 
2mLo^ 



V's(r,v) 



, cxp < — -X ■ 

(27r)2VdetM I 2 



(Al) 



covariance matrix 



M 



_ ( {YV^)n (rv'^>o 
(vr'^)o (vv'^)o 



alone. We calculate M using Chandrasekhar's theo- 
rem [40], 



M 



™^ ./o 



At 



G™(i)G™T(t) G™(i)G™T(i) 
G™(t)G™T(i) G™(t)G™T(t) 



With 1 + A{mioY — Akm, we obtain the stationary cor- 
relations 



1 







(rv^>o 



1 

^ / -7 
2/c I 7 



l+4fcm -2 • 

7 



(A2) 
(A3) 
(A4) 



The distribution function projected into the configura- 
tion space is obtained through integrating out the veloc- 
ities in the stationary distribution (Al), 



= exp s — -r ■ (rr )n r 

2^v/det(rrT)o ^12^ ^" 



In the overdamped limit one obtains 



(rr^)o 



>o 1 / 1 -|-_272 7 
k[ 1 1 



is Gaussian and therefore determined by the symmetric with 7 — ^ /2k. 
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